RL-TR-96-240 
Final  Technical  Report 

February  1997 


THEORETICAL  FOUNDATION  OF  HIGH- 
FREQUENCY  MICROWAVE  MODULATION 
OF  LASER  DIODES  WITH  SHORT 
EXTERNAL  CAVITY 


Rochester  Institute  of  Technology 
Dr.  Guifang  Li 


DTTc 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


19910416  006 


Rome  Laboratory 
Air  Force  Materiel  Command 
Rome,  New  York 


This  report  has  been  reviewed  by  the  Rome  Laboratory  Public  Affairs  Office 
(PA)  and  is  releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RL-TR-96-240  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


‘inn 


/ 


STEVEN  JOHNS 


Project  Engineer 


FOR  THE  COMMANDER; 


DONALD  W.  HANSON,  Director 
Surveillance  &  Photonics  Directorate 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  Rome  Laboratory 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  please 
notify  RL/OCPC,  25  Electronic  Pky,  Rome,  NY  13441-4514.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a  specific 
document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


1 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  compieting  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson 
Havk  Hinhwav  Suite  1204.  Arlinoton.  VA  22202-4302,  and  to  the  Office  of  Manaqement  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1  .  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

February  1997 

3.  REPORT  TYPE  AND  DATES  COVERED 

FINAL,  May  95  -  Dec  95 

4.  TITLE  AND  SUBTITLE 

THEORETICAL  FOUNDATION  OF  HIGH-FREQUENCY  MICROWAVE 
MODULATION  OF  LASER  DIODES  WITH  SHORT  EXTERNAL  CAVITY 

5.  FUNDING  NUMBERS 

C  -  F30602-95-C-0091 

PE  -  62702F 

PR  -  4600 

TA  -  P2 

WU  -  PU 

6.  AUTHOR(S) 

Dr.  Guifang  Li 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES| 

Rochester  Institute  of  Technology 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

75  Highpower  Rd 


Rochester  NY  14623 


N/A 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Rome  Laboratory/OCPC 
25  Electronic  Pky 
Rome  NY  13441-4515 


1 1 .  SUPPLEMENTARY  NOTES 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 

RL-TR-96-240 


Rome  Laboratory  Project  Engineer:  Steven  T.  Johns,  OCPC,  (315)  330-4456 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  Public  Release;  Distribution  Unlimited 


13.  ABSTRACT  (Maximum  200  words) 

In  simulations,  short  external  cavity  laser  diodes  have  shown  to  produce  high  frequency  oscillations 
(>20Hz).  This  paper  presents  theoretical  and  simulated  dynamics  for  such  types  of  laser  systems.  The 
high  frequency  oscillations  are  a  result  of  self  twin-mode  locking.  It  will  be  shown  the  critical  feedback 
parameters  and  the  frequencies  of  self  twin-mode  locking  oscillations  in  these  regimes  of  high  frequency 
microwave  oscillations  are  possible. 


14.  SUBJECT  TERMS 

high  frequency  microwave  modulation,  short-external  cavity, 
semiconductor  laser  diodes 

15.  NUMBER  OF  PAGES 

52 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UNLIMITED 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 


Prescribed  by  ANSI  Std.  Z39-1 8 
298-102 


Table  of  Contents 


Abstract  3 

1.0  Introduction  4 

2.0  Theory  6 

2.1  Rate  Equation  6 

2.2  Steady-State  Analysis  7 

2.3  Compound-Cavity  Modes  9 

2.4  Linear  Stability  Analysis  15 

3.0  Numerical  Simulation  16 

3.1  F abry-Perot  Laser  with  an  External  Cavity  1 6 

3.2  VCSEL  with  an  External  Cavity  23 

4.0  Conclusion  28 

Appendix  A  29 

Appendix  B  31 

Appendix  C  34 

References  43 


1/2 


Abstract 


The  dynamic  properties  of  semiconductor  laser  diodes  are  significantly  affected  by 
external  optical  feedback.  In  this  paper  we  present  theoretical  and  simulated  dynamics  for  a  laser 
diode  with  a  short  external  cavity.  We  assume  that  only  the  emitted  photons  are  fedback  only 
once  by  the  external  cavity.  In  the  unstable  regime  of  short  external  cavity  laser  diode  operation 
high  frequency  oscillations  (>20  GHz)  exist  as  a  result  of  self  twin-mode  locking.  This  instability 
corresponds  to  equalization  of  the  cavity  loss  for  the  main  mode  (the  lasing  frequency)  and 
sidemode.  We  will  show  the  critical  feedback  parameters  and  the  frequencies  of  self  twin-mode 
locking  oscillations  in  these  regimes  of  high  frequency  microwave  oscillations.  As  will  be 
shown,  both  the  Fabry-Perot  and  vertical  cavity  surface  emitting  lasers  exhibit  such  microwave 
oscillations.  The  frequency-modulated  microwave  oscillation  will  be  useful  for  all-optical  high¬ 
speed  transmitters. 
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1.0  Introduction 


Semiconductor  laser  diodes  exposed  to  optical  feedback  have  been  found  to  improve  laser 
frequency  stability  and  to  narrow  laser  linewidth.  Optical  feedback,  however,  can  also 
significantly  change  LD  (laser  diode)  dynamics. 

It  also  has  been  theoretically  and  experimentally  demonstrated  that  the  effects  optical 
feedback  vary  with  the  external  cavity  length  and  with  power  feedback  ratio.  For  long  (  ~lcm 
and  more)  external  cavity  LD,  these  changes  lead  to  chaotic  self-modulation  of  laser  light  with 
dramatic  increase  in  laser  spectrum  width  and  noise,  instead  of  line  narrowing  (  coherence- 
collapse  “  effect  [1]).  It  occurs  as  a  result  of  sustained  relaxation  oscillations  of  laser  output, 
when  feedback  coefficient  exceeds  some  very  low  critical  value  about  -40  dB  and  extends  over 
several  decades.  This  “coherence-collapse”  regime  was  shown  to  be  the  result  of  a  combined 
effect  of  the  LD  relaxation  resonance  and  EC  resonances  [2]  .  For  the  short  external  cavity  LD 
(SEC  LD),  the  EC  mode  spacing  is  much  larger  than  the  LD  relaxation  frequency  [4].  The  larger 
mode  spacing  lend  to  avoid  coherence-collapse. 

In  the  unstable  regime,  a  SEC  LD  exhibits  high-frequency  instability  associated  with 
locking  of  compound  cavity  modes  for  relative  high  values  of  coupling  coefficients.  This  effect 
enables  to  achieve  frequencies  of  optical  current-controlled  oscillators  in  excess  of  20  GHz. 

We  can  use  SEC  LD  to  generate  tunable  microwave  oscillators  with  frequency 
modulation  capacities.  The  frequency  can  be  modulated  by  several  ways.  From  our  simulation 
result,  we  can  modulate  the  frequency  by  changing  optical  feedback  and  the  external  cavity 
length. 


In  the  following,  we  focus  first  on  the  theory  of  SEC  LD  including  rate  equation,  steady 
state  analysis  and  linear  stability  analysis.  We  introduce  the  Lang  and  Kobayashi  rate  equations 
for  the  SEC  LD  [3].  The  steady  state  analysis  shows  that  for  each  state  of  lasing  the  “steady- 
state”  condition  of  instability  occurs  when  the  cavity  losses  for  the  main  and  side  mode  are  equal 
[4].  Also,  the  linear  stability  analysis  can  predict  the  dynamic  instability  of  laser  operation  caused 
by  the  twin-mode  competition. 

Furthermore,  we  use  computer  simulation  to  prove  the  instability.  First,  Fabry-Perot  laser 
diodes  with  200  urn  cavity  length  is  used  for  computer  simulation.  Second,  we  use  vertical-cavity 
surface  emitting  lasers  (VCSEL)  with  1  um  cavity  length  for  simulation.  The  structure  of 
VCSEL  is  shown  in  Fig.  1  [5]. 

Compared  with  long  cavity  length,  the  advantage  of  the  VCSEL  is,  due  to  its  short  cavity 
length,  the  longitudinal  cavity  mode  spacing  (being  equal  to  c/2d  ,  d  is  the  cavity  length)  is  very 
large.  In  fact,  it  is  possible  to  have  a  mode  spacing  larger  than  the  gain  bandwidth  of  the  active 
material.  Thus,  it  can  be  inherently  single  longitudinal  mode. 
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Finally,  our  results  are  summerized. 
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Fig.  1  Schematic  cross-section  view  of  VCSEL  structure  considered  here. 
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2.0  Theory 

2.1  Rate  Equation 


Laser  Diode  External  Cavity 


R1  R2  R3 


1 ^ — 

L  Lext 


Fig.  2  Schematic  of  a  laser  diode  with  an 
external  cavity 

The  schematic  of  a  semiconductor  laser  diode  with  an  external  cavity  is  shown  in  Fig. 2  . 
R1  and  R2  are  the  reflectivity  of  the  laser  facets  and  R3  is  the  reflectivity  of  the  external  mirror. 
L  and  Lext  are  the  lengths  of  the  laser  diode  cavity  and  the  EC;  fext  is  the  feedback.  We  will 
consider  xco^j  «1  ,  where  r  =  2 1^,^,  /c  is  the  photon  round-trip  time  in  EC  and  co^  is  the 
angular  frequency  of  the  relaxation  oscillation  . 


Here  we  use  Lang  and  Kobayashi  rate  equations  for  the  SEC  LD.  There  are  several 
assumptions  for  these  rate  equations.  First,  the  solitary  laser  is  assumed  to  be  single  mode,  which 
implies  r»  (tq  =  'ILn^  /  c  ,  is  the  photon  round-trip  time  in  the  LD-cavity.  is  the  group 
refraction  index  of  the  LD.)  ;  second,  «  1  which  means  only  one  round-trip  of  the  emitted 
photons  is  significant  inside  the  EC.  Then  we  arrive  at  the  following  two  rate  equations  : 


dm 

dt 


E(r)  +  — (g — ^)£'(0 +  /cE(r  -  •r)exp(-iQr) 
2 


(1) 


dt 


(2) 


where  g  =  g{N){\-  k^P)  is  the  nonlinear  gain  where  the  linear  gain  is  given 
by  W)  =  gjv  ( iV  -  ) ;  P  is  the  output  power. 

is  the  resonance  optical  frequency  of  the  LD  without  EC. 
a  is  the  linewidth  enhancement  factor;  g,^  is  the  differential  gain. 
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Tp  is  the  photon  lifetime  in  the  LD  cavity. 

is  the  carrier  lifetime  due  to  spontaneous  and  nonradiative  recombination. 

iS  =  |£(t)|^  is  the  photon  number  in  the  LD  cavity,  (output  power  P~S) 

J  is  the  carrier  injection  rate. 

The  feedback  parameter  is  given  by 


where  parameter  C^  which  is  a  function  of  the  facet  and  external  mirror  reflectivity  and  R.^ 


V  ^2 

is  a  measure  of  the  coupling  strength  between  the  two  cavities. 

The  relation  between  output  power  (P)  and  photon  number  (S)  is  discussed  more  in  detail 
in  the  appendix  A  [6]. 

2.2  Steady-state  analysis 


Stationary  lasing  conditions  can  be  obtained  from  (1),  by  setting  E  to  be  constant.  The 
real  and  the  imaginary  parts  of  the  equation  ,when  they  are  separated,  can  be  written  as 


g — ^  =  -2/ccosQr 


a 


-Q  +  —  Ag- A:sinQr=  0  . 


(3) 

(4) 


In  the  absence  of  feedback,  the  threshold  gain  g,^  is  determined  by 

R,h=—  ■  (5) 

Inserting  (5)  into  (3)  leads  to  a  threshold  gain  shift 

=  =-2/rcosQr.  (6) 

To  find  possible  frequencies  of  lasing  ,  insert  (6)  into  (4) ;  that  gives 
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-  Q  =  /cacosQr+  /csinQr 


(7) 


or 

(Q^  -Q)r=  A:T(acosQ7+sinQr)  =  KT^ll  +  sin(Qr+  ff) 

{Q.^-0.)t=  Csin{O.T+  ff),  (8) 

where 

0  =  arctan  a  (9) 

and  the  feedback  coefficient  C  is  given  by 


C  =  KT^l\  +  a^  .  (10) 

The  number  of  external-cavity  modes  is  determined  by  the  parameter  C.  For  C<1  only 
one  mode  exists,  while  for  C>1  there  exists  two  or  possibly  more  modes.  For  C»1  the  total 
number  of  modes  is  approximately  given  by  C/{7r+\)  (Fig. 3  shows  the  solutions  when 
increasing  kt).  Furthermore  for  C>1,  only  half  of  the  solutions,  corresponding  to  dCll  dQ^  >  0 
are  stable  with  respect  to  slow  fluctuations,  representing,  thus  ,  possible  frequency  of  stationary 
lasing.  Qq  denotes  the  possible  frequency  of  lasing.  Inserting  Qp  into  Q  of  equation  (6)  and  (8), 
gives: 


Fig.  3  The  graphical  representation  of  equation  (11)  for  a=4.  dashdot  line  for  /cn=0.5  , 
solid  line  for  and  dash  line  for  /cr==l 6  .  intersection  between  the  three  lines 

and  sinusoidal  function  is  theoptical  frequencyof  the  LD  without  EC. 
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Vl  +  cc^  sin(Qo  r+  arctan  a) 


(11) 


K 

^^  =  -cosQor. 
2k 


(12) 


The  graphical  representation  of  equation  (11)  is  shown  in  Fig.  3.  The  solutions  to 
equation  (11)  are  the  intersections  between  the  line  and  the  sinusoidal  function  [7].  Also  when 
Kr<\  only  one  solution  exists  (the  intersection  between  the  dashdot  line  and  the  solution  of  the 
LHS  of  equation  (11))  ,  increasing  kt  becomes  more  solutions  (  shown  by  the  intersections 
between  the  other  lines  and  the  sinusoidal  function) . 


2.3  Compound-Cavity  modes 

Each  lasing  frequency  Qq  has  its  own  system  of  compound-cavity  modes  Q„.  The 
lasing  frequency  Qq  and  the  compound  cavity  modes  should  have  the  same 


normalizid  optical  frequency  shift 


Fig.4  normalized  cavity  loss  shift  versus  normalized  optical  frequency  shift,  the  ellipse 
correspond  to  possible  lasing  frequency.  The  circle  shows  CC  modes  when  the  lasing 
frequency  is  equal  to  the  optical  frequency  of  the  LD  without  EC. 

Ag,/j(  threshold  gain  shift  due  to  feedback  ).  That  gives 

=  ^S,h  (^o)  =  -iKcosa^T.  (13) 
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By  inserting  (13)  and  Q„  into  (7),  we  obtain 


-Q„  =  ^sinQ„T-yAg(Qo) .  •  (14) 

These  compound-cavity  modes  Q„  can  be  found  by  (14)  times  r  and  subtracting  the  result  from 
(11)  times  KT ,  which  yields 

— ^^  =  sinQ„r-sin^o  .  (15) 

K 

The  graphical  representation  of  (15)  is  shown  in  Fig  .  5.  The  solutions  to  (15)  are  the 
intersections  between  the  line  and  the  sinusoidal  function.  As  we  can  see  in  Fig.  5,  increasing  kt 
cause  more  and  rnore  solutions  between  the  LHS  and  RHS  of  equation  (15). 

Let  us  consider  the  evolution  of  the  CC  modes  with  increasing  kt  in  the  plane  of 
parameters  normalized  optical  frequency  shift  (Qo-Q^)/a:  --  normalized  cavity  loss  shift 
^Sioss  /2k=  -  cosOqT.  In  this  plane,  see  Fig.  4  ,  all  solutions  of  (8)  belong  to  an  ellipse  [8].  For 
Kr»  1  only  those  solutions  of  (11)  that  lie  on  the  lower  branch  of  the  ellipse  are  unstable 
corresponding  to  JQq  /  dQ.^  <  0  [2]. 

The  solutions  of  (15)  lie  on  a  circle  with  radius  1  in  the  . same  coordinates  as  Fig.  4. 
Equation  (8)  AQr=  Qr-QoT= -Csin(Qr+arctana)  and  a  small-signal  analysis  yield  the 
result  for  the  linewidth  of  narrowing  by  external  feedback  [12]: 


Av=  Avo[l  +  Ccos(Qr+arctana)]  (16) 

where  A  Vg  denotes  the  linewidth  of  the  solitary  laser  diode. 

The  lowest  spectral  linewidth  is  obtained  for  =  Qr  =  -arctan  a,  yielding 


Av: 


Av„ 


AVn 


(17) 


(1  +  C)^  (l+KrVTW)^  ■ 

The  mode  with  the  lowest  threshold  gain  to  start  lasing  is  Ag,/,  = -2KCOsQr= -2k 
when  cos  Q  r  =  1  and  thus 


Av: 


Av„ 


(l  +  srr)^ 


for  (j)^^,=27tm  ,  7w  =  0,1,2,- 


(18) 
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Fig.5  The  graphical  representation  of  equation  (15).  Dash-dotted  line  for  Kr=Q.5,  solid 

line  for  kv=A  and  dash  line  for  a:z=16. 

One  would  more  expect  that  the  mode  with  the  lowest  threshold  gain  would  be  the  lasing 
one.  Surprisingly,  it  turns  out  that  the  measured  spectral  line  width  for  laser  diode  with  optical 
feedback  coefficients  C»1  is  closer  to  equation  (16)  than  to  (18),  indicating  that  the  external 
cavity  mode  with  the  narrowest  linewidth  is  more  stable  than  the  external  cavity  mode  with  the 
lowest  threshold  gain. 

From  the  measured  conclusion,  we  assume  QQT=-arctana  corresponds  to  laser 
oscillation  with  minimal  linewidth.  In  Fig.  4  the  circle  is  shown  for  the  case  of  Qq  =  Q,  .  It 
means  that  the  lasing  frequency  Qq  is  equal  to  the  resonance  optical  frequency  of  the  LD 
without  EC.  The  intersections  between  the  ellipse  and  the  circle  are  four  points.  One  of  the  points 

marked  A  in  Fig.  4  is  when  Qq  =Q,  and  its  coordinates  are  [(Qq  -  Qj//c,-cos{^oj.  Another 
point  B  has  coordinates  ,  [(Qi  -  QJ//c,-cos^o]  where 

Qj  =  Qo +2/csin^o  .  (19) 

The  points  A  and  B  have  the  same  optical  loss  =  -2/ccosQor  as  the  lasing  mode 
Qg ,  they  can  be  oscillating  modes  at  the  same  time. 

For  KT<  1  (weak  feedback),  Qg  =  Q„  is  the  only  solution  to  equation  (15).  And  it  cannot 
be  the  lasing  one  because  the  optical  loss  is  higher  than  the  lasing  one  (point  A).  Increasing  kt,  it 
leads  to  equalization  of  the  cavity  losses  for  the  main  mode  (  point  A)  and  the  sidemode  when 
the  last  one  reaches  point  B.  Point  B  (the  sidemode)  has  the  same  optical  loss  as  point  A  (the 
main  mode)  and  could  be  the  lasing  one  too.  However  point  B  is  lying  on  the  upper  unstable 
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branch  of  the  ellipse,  so  it  cannot  be  the  lasing  one  for  the  solitary  LD  without  EC.  For  the  LD 
with  EC,  reaching  point  B  by  the  CC  sidemode  can  be  considered  as  the  transition  from  stable 
CW  operation  at  frequency  Qq  to  instability.  The  critical  value  of  kt  to  reach  point  B  can  be 
found  by  the  following  steps.  From  the  x-axis  of  Fig.  4  at  point  B, 

cosQiT=  cos^zIq  s  cosQoT  .  (20) 

Q,  also  meets  one  of  the  solutions  of  equation  (15).  By  inserting  Qj  into  of  equation  (15), 
we  get 


(Qq  - Q,)z-  KTisinCl^r -sin (po) 


(21) 


Inserting  equation  (19)  Q;  =  +2Ksmpf^  into  Q,  on  the  LHS  of  equation  (21)  yield 

sinQ; r= -sin^o  .  (22) 

From  equations  (21)  and  (22),  we  have 


Q,r=  Qor-2;T-2^Q  .  (23) 

After  multiplying  equation  (18)  by  t  and  subtracting  (22)  into  it,  one  obtains  the  “steady-state” 
condition  of  instability  of  the  form 


KT- 


-  sin  ^0 


Equation  (24)  also  satisfies  equation  (11). 


(24) 


Equation  (24)  defines  the  critical  value  of  feedback  under  two  conditions: 
=  -arctana  and  Qp  =  .  We  could  associate  each  state  of  lasing  with  its  own  number  N. 

The  point  A,  one  of  lasing  state  nearest  is  set  to  be  N=0.  We  denote  Qp''^,  the  lasing 
frequency  for  the  Nth  state,  and  ^zip^ ,  the  detuning  parameter  for  the  Nth  state.  Then  we  can  find 
the  next  relation  for  the  phase  shift 


Qp""  =Q,r-2;iiV-^,  +  {zip"'  iV  =  0,±l,±2, . .  (25) 

<f)  is  the  detuning  parameter  between  {-n,7r)  ,  and  £^.r=2nn  +  . 

Combining  (24)  and  (25)  in  (7),  we  obtain 
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“  7t{2N  +  \)^(l), 


(26) 


This  is  the  equation  for  the  critical  value  of  phase  detuning  . 

To  find  the  critical  feedback  coefficient  value  to  reach  the  “steady-state”  condition  (24)  of 
instability  for  each  state  of  lasing,  we  combine 


1-^  r  fy 


for  Fabry-Perot  laser  and  (23).  That  gives 


Kt  = 


\-R  T 


■iR  I'd 


-  sin  ^0 


then, 


J  ext 


R 


+  f  T 


(l-Ry 


(27) 


The  relation  between  main  mode  and  sidemode  (CC  mode)  of  the  Nth  lasing  state  can  be 
found  by  using  equation  (23) , 

i.e.. 

Then, 
where 


7t  X 


Qj^  =  -  InVff 


(28) 


= 


1-h 


^0 


y 

X 


(29) 


V  ^7 

can  be  interpreted  as  the  beating  frequency  between  two  CC  modes  having  equal  optical  loss. 

In  equation  (24)  kx  should  be  greater  than  zero,  so  has  values  between  (-/r,0). 
That  is,  the  detuning  parameter  is  negative.  This  condition  translates  to 


^0 


N 


1+^<1 

n 


which  gives 
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(30) 


1 

^N<- 
X 

from  equation  (29).  So  the  mode  beating  frequency  is  always  less  than  EC-mode  spacing  1  /  r. 

The  steady-state  analysis  of  the  short  EC  LD  will  help  us  to  find  the  critical  value  of  the 
feedback  parameter  leading  to  possible  instabilities  in  laser  operation.  Now,  we  will  prove  the 
instability  of  laser  operation  by  using  the  dynamic  analysis  including  the  linear  stability  analysis 
and  computer  simulation. 

2.4  Linear  Stability  analysis 

The  static  solution  is  called  stable,  if  small  perturbations  from  the  static  solution  decay  to 
zero  as  time  goes- toward  infinity.  Such  a  stable  static  solution  can  be  observed  experimentally. 
To  check  the  stability  of  the  static  solution,  we  perform  the  linear  stability  analysis  for  a  given 
static  solution.  Let  5E,5^  ,and  5N  denote  the  deviation  from  the  stationary  solution  .  The 
criterion  for  stable  operation  is  that  the  determinant  D(s)  has  no  zeros  in  the  right  half  of  the  S- 
plane  [7],  [9].  The  determinant  is  obtained  as 

D{s')  -  ~s^  +  s^[y  +  2 kb{t) cos <f)\ 

-s^cOf^  +  {t)  +  2y^kb{t)cos^  (31) 

+  }^co,^^e{T){cos^-  a  sin  (fi)  +  Y  +  crXeXp^(^)sin^zl] 

with  c(  r)  =  1  -  and  <j)  =  Q.T.  Here  =  1  /  is  the  carrier  perturbation  decay  rate, 

y-Ye’^Yp  is  the  decay  rate  of  relaxation  oscillations,  with  Yp  =  g^pP  being  the  contribution 
of  the  nonlinear  gain.  The  angular  frequency  of  relaxation  oscillations  (Op  is  given  by 
o>R  Ye  Yp  ■  Assuming  that 

^/»YeYp  ,  G>R»J<Ye 

then  D(S)  can  be  approximated  by 

D{s)  =  -s^  +5^[y  +  27CB(r)cos<zi] 

- s^o) p  +  (32) 

+  j<^cOp^e{T){cos(^- asin  ^)j 
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For  s  =  ja>,  we  can  separate  the  equation  (32)  into  the  real  and  imaginary  parts.  The  real  part  of 
the  equation  is  related  to  the  damping  constant  of  the  semiconductor  laser  rate  equations,  and  the 
imaginary  part  is  related  to  the  small  signal  resonance  frequency. 

The  real  and  the  imaginary  parts  give  the  following  two  conditions  for  having  stable 
solutions  located  in  the  left  half  of  the  complex  S  plane  (The  detail  is  shown  in  Appendix  B) : 

(33) 

(34) 

where  co^  =  co^  /  y  is  the  damping  frequency  of  the  relaxation  oscillations  of  LD  without  EC, 
(f)^  =  [Qq  r]  mod (2  n)  is  the  optical  phase  shift  in  EC,  or  detuning  parameter. 

If  we  neglect  the  first  term  on  the  right  hand  side  (RHS)  of  the  equations  (33)  and  (34), 
the  remaining  terms  H  and  y  describe  the  relaxation  oscillation  -  induced  instability  which 
causes  coherence  collapse  for  long  EC  LD  and  very  weak  feedback.  The  first  terms  in  the  RHS 
of  equation  (33)  and  (34),  which  are  proportional  to  and  k,  respectively,  are  thus  of 
important  in  the  case  of  sufficiently  high  levels  of  feedback  and  reflect  the  influence  of  the 
above-described  CC-mode  competition  on  stability  of  the  short  EC  LD. 

In  the  case  of  sufficiently  high  levels  of  feedback,  we  neglect  the  second  term  on  the  RHS 
of  both  (33)  and  (34)  and  assume  2o}^  /  » \,a.  We  then  obtain  fi’om  (33)  and  (34) , 


n  T  A  ■  A 

(1-2 — ^)cos^o  -orsm^o  <2 — j  sin  a)z+ 

2  2  A  2  •  2  /  /■  OyT. 

CO  -cOj^  =Ak  sm  (— )  +  (S)ycot(— ) 


KO)j  1  -  cos  COT 


A  ^  • 

COS  - sm  COT 

CO 

CO  =  2/csin(-^) 


(35) 

(36) 


When  the  frequency  of  self-oscillations  col 2k  is  equal  to  the  CC-mode  beat  frequency  from 
equation  ( 1 9),  we  get 


CO  1 

—  =  v„=(l  +  ^)-- 
2k  k  t 


(37) 


which  yields 


cot=2k^-2^^  (38) 

By  inserting  (38)  into  (35)  or  (36)  ,we  can  get  the  same  expressions  as  (14) : 
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-sin 

Thus,  the  linear  stability  analysis  of  the  short  EC  LD  not  only  gives  the  critical  value  of  k  but 
also  predicts  the  dynamic  instability  of  laser  operation  caused  by  the  competition  of  the  two  CC 
modes  which  have  equal  optical  losses. 

Equations  (31)  and  (32)  give  the  boimdary  of  this  stability  region.  But  they  cannot  give 
any  information  above  the  stability  boundary.  To  get  the  information  above  the  boundary  of  the 
stability  region,  we  proceed  to  use  computer  simulations  of  the  nonlinear  system. 


3.0  Numerical  Simulation 


The  rate  equations  for  the  EC  LD  are  given  in  equations  (1)  and  (2).  The  absolute  squared 
value  of  the  electric  field  amplitude  corresponds  to  the  number  of  photons  in  the  laser  cavity. 
Equations  (1)  and  (2)  can  be  cast  into  rate  equations  for  the  photon  number  S,  the  phase  (f)  ,  and 
the  carrier  number  N. 


dSjt) 

dt 

d(t>{t) 

dt 


=  (g-—)S(t)  +  2K^Sit).JS(t-T)  cos(Q,r+  ^(t)  -  (f>{t  -  r)) 
= -^^iv  (^  -  -  ^^^^^^^sin(Q,T+ (/5(0  -  <z^(t  -  r)) 


(39) 

(40) 


dN{t) 

dt 


-sp 


(41) 


In  numerical  simulations,  the  system  of  rate  equations  system  (39)-(40)  is  solved  by  using 
a  fifth-order  Runge-Kutta  method  [10],[1 1].  The  calculations  used  the  laser  parameters  given  in 
Table  I  .The  computer  program  is  written  in  C  (shown  in  Appendix  C).  We  specify  the  initial 
conditions  S(0)=10,  ^zi(0)=0  and  N(0)=0.  In  the  time  interval  0<r<  t,  the  data  of  photon 
number  S  and  phase  (f)  are  stored  in  memory  and  they  are  reused  after  the  round-trip  time  r  of 
the  EC  in  order  to  account  for  the  optical  feedback. 


Next,  we  will  show  results  for  the  two  different  kinds  of  SEC  LD:  A)  Fabry-Perot  laser 
with  an  external  cavity;  B)  VCSEL  with  an  external  cavity 


3.1  Fabrv-Perot  Laser  with  an  external  cavity 

The  first  example  is  a  Fabry-Perot  laser  with  an  external  cavity  as  shown  in  Fig.2,  with 
cavity  length  L=200um,  facet  reflectivity  R=0.32  (  corresponding  to  KL=2.3  for  an  AR-coated 
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DFB  laser  with  di  nil  phase  shift.),  output  power  P=5mW/per  facet.  Other  parameters  are  given 
in  Table  I  [12].  We  choose  EC  cavity  length  such  that  round-trip  time  z  varies  between  25  and 
40  ps  (the  external  cavity  between  3.75  and  6  mm).  Within  this  range,  r  is  less  than  the  critical 
value  for  coherence-collapse  regime  (chaotic  self-modulation  of  laser  light  with  dramatic 
broadening  of  laser  spectrum). 

TABLE  I 

List  of  Parameters  for  Fabry-Perot  laser 


Parameter 

Symbol 

Value 

Units 

scattering  loss  in  the  active  region 

3.6x10^ 

m-' 

carrier  number  at  transparency 

Ns 

2.0x10' 

output  power 

P 

2.906x10''  X  5 

mW  /  facet 

index  of  active  region 

4.5 

differential  gain 

Sn 

2.2x10' 

s-' 

linewidth  enhancement  factor 

a 

4 

laser  wavelength 

X 

1.3 

um 

spontaneous  emission  lifetime 

2.0x10-' 

s 

photon  lifetime 

1.29072x10'" 

s 

gain  saturation  coefficient 

K 

4.8 

w-' 

Figs.  6  and  7  illustrate  two  characteristics  of  SEC  LD  instability.  In  Fig.  6  we  report  the 
simulated  and  analytical  dependence  of  the  oscillation  according  to  the  feedback  parameter  (/„, ) 
on  the  length  of  external  cavity  (Text).  The  dependence  of  the  frequency  of  these  oscillations  on 
the  length  of  external  cavity  (Text)  is  shown  in  Fig.  7  .  The  two  dash  lines  in  Fig.  6  and  Fig.  7 
represent  two  groups  of  strongest  high-frequency  oscillations  (the  modulation  depth  reaches 
unity)  according  to  numerical  simulation.  Data  marked  with  asterisks  are  the  simulation  data  of 
high  frequency  oscillations.  The  lower  dash  line  belongs  to  the  first  group  of  high  frequency 
oscillation  according  to  .The  second  dash  line  belongs  to  the  second  group  of  high 
frequency  oscillation  according  to  . 

As  discussed  the  results  in  Fig.  6  and  Fig.  7  are  obtained  under  two  assumptions.  First, 
=  -arctancr,  the  lowest  spectral  linewidth,  is  more  stable  than  the  lowest  threshold  gain  as 
mentioned  before.  Second,  means  the  lasing  frequency  Qq  is  equal  to  ,  the 

resonance  optical  frequency  of  the  LD  without  EC  .  We  can  see  from  the  high  feedback 
parameter  (/^^,)  that  the  laser  operation  will  lose  stability  and  lead  to  high  frequency 
oscillation. 

This  type  of  instability  is  the  result  of  the  beat  frequency  between  two  CC  modes  having 
equal  optical  loss.  As  mentioned  in  the  previous  section,  by  increasing  /^,  it  leads  to 
equalization  of  the  cavity  losses  for  the  main  mode  and  the  first  left  sidemode  and  those  two 
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modes  produces  a  self  twin-mode  locking  generating  high  frequency  oscillations.  In.  Fig  6  the 
two  solid  lines  present  the  critical  feedback  parameters  from  equation  (27)  (the  lower  solid  line 
shows  N=0  for  and  the  other  shows  N=l),  corresponding  to  the  main  mode  and  the 

sidemode  having  equal  optical  loss.  The  two  solid  lines  as  shown  in  Fig.7  illustrate  the  results 
from  equation  (29)  (the  lower  solid  line  shows  N=0  for  and  the  other  shows  N=l).  These 
results  denote  the  beat  frequencies  between  two  CC  modes  having  equal  optical  loss.  From  Fig.  6 
and  Fig.  7  linear  stability  analysis  and  computer  simulated  agree  very  well. 


Fig.  6  A  SEC  LD  operation  (Lext,fext).  The  lower  and  upper  solid  lines  present 
theoretical  critical  value  at  and  group  .  The  lower  and  upper  dash  lines  show 
the  strongest  oscillation  at  and  group.  Curved  marked  with  represent  the 
simulation  data.. 

In  what  follows,  we  explain  the  data  within  the  onset  and  offset  of  high  frequency 
oscillations  of  the  group  for  an  air-equivalent  EC  length  Lext=5.25mm  (t=35  ps).  For 
=0.204  the  high  frequency  oscillation  occurs  at  a  frequency  around  20.9  GHz.  But  the  oscillation 
is  not  strong,  the  amplitude  of  the  oscillation  will  become  weaker  and  it  will  disappear  after 
t>52ns.  Increasing  to  0.205  the  stable  oscillation  occurs.  But  the  amplitude  of  the  oscillation 
is  small.  That  means  the  intensity  of  the  sidemode  is  much  less  than  the  main  mode.  Further 
increase  in  /^,  leads  larger  oscillation  amplitude  (the  stronger  intensity  of  the  sidemode)  and  the 
bigger  oscillation  frequency.  When  reaches  0.206  the  strongest  intensity  shows  up  as  shown 
in  Fig.  8  with  the  frequency  f=21.3  GHz.  Further  increase  in  /^,  (>0.206)  destroys  the  twin¬ 
mode  locking.  It  leads  the  high  frequency  oscillation  to  a  new  single-mode  state  with  higher 
output  power  (6.56  mW).  The  new  optical  frequency  does  not  the  same  with  any  of  the  “old”  CC 
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modes  in  Fig.  8  .  The  new  one  is  one  of  the  stable  solutions  (which  lie  on  the  lower  branch  of 
Fig.  4)  of  equation  (8). 


Fig.  7  A  SEC  LD  operation  (Lext,f).  The  lower  and  upper  solid  lines  present  the 
theoretical  critical  value  at  and  group  .  The  lower  and  upper  dash  lines 

show  the  beat  frequencies  of  the  and  according  to  the  strongest  oscillation 
in  Fig.  11.  Curved  marked  with  **’  represent  the  simulation  data. 

From  the  theoretical  and  simulated  results,  we  conclude  that:  for  the  high  feedback 
parameter  ,  increasing  causes  the  single-mode  optical  frequency  to  become  unstable 
above  a  critical  .  It  leads  to  high-frequency  oscillation  corresponding  to  self  twin-mode 
locking.  Increasing  f^,  destroys  mode  locking.  If  feedback  /^,  is  further  increased,  above 
process  will  repeat. 

Thus,  we  can  see  that  a  SEC  LD  exhibits  several  regimes  of  self  twin-mode  locking  high 
frequency  oscillations  in  some  ranges  of  feedback  parameters. 

The  simulation  we  have  done  corresponds  to  an  output  power  P=5mw/per  facet.  It  agree 
with  the  results  in  [4].  The  critical  feedback  parameters  required  for  high  frequency  oscillation 
also  change  with  the  change  of  the  injection  number  according  to  the  simulated  results. 
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Fig.  8  The  strongest  oscillation  for  the  Lext=5.25mm  and  the  group.  The 
feedback  parameter  is  0.206.  The  oscillation  frequency  is  equal  to  21.3  GHz. 

In  what  follows,  we  set  the  output  power  P=2  mW/per  facet.  Other  parameters  are  the 
same  as  the  previous  example  as  shown  in  Table  I.  Here,  we  also  choose  the  EC  length 
=5.25  mm.  We  found  that  the  feedback  parameters  for  the  strongest  modulation  index 
occurs  when  /„,=  0.188  compared  to  0.206  for  P=5  mW/per  facet.  Further  increasing 
destroys  the  twin  mode  locking.  It  leads  the  laser  operation  on  a  new  single-mode  state  with 
higher  output  power  (3.06  mW).  The  strongest  oscillation  frequency  we  get  in  the  group  /„/  is 
20.7  GHz.  This  frequency  is  not  the  same  with  the  frequency  21.3  GHz  in  the  previous  example 
of  P=5  mW/per  facet. 

Table  II 

List  of  frequency,  average  power  and  oscillation  amplitude 
on  the  change  of  the  injection  current 


Injection  current  (mA) 

frequency  (GHz) 

Average  Power  (mW) 

Osc.  Amp.  (mW) 

47.87 

0 

6.279 

0 

48.23 

21.3 

4.925 

4.775 

48.39 

21.2 

4.95 

4.95 

48.55 

21.1 

5.06 

4.89 

48.71 

21.05 

5.05 

4.55 

48.84 

21.03 

5.12 

4.0 

48.85 

21.02 

5.123 

3.925 

48.865 

21.01 

5.125 

3.825 

48.87 

21 

5.163 

0.123 

49.03 

20.9 

5.228 

0.007 

49.19 

20.7 

5.2875 

0.0005 
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Next,  we  will  show  the  dependence  of  the  oscillation  frequency  (f),  average  power  (P) 
and  oscillation  amplitude  on  the  injection  current  (I=J/e).  The  external  cavity  length  is  fixed  at 
Lext=5.25  mm  and  the  feedback  parameter  is  fixed  at  =0.206.  The  other  parameters  are 
shown  in  Table  I.  In  Table  II  and  Fig.  9  (a),(b),(c)  ,  we  can  see  the  lead  to  higher  oscillation 
frequency  variations  of  these  quantities  for  injection  current  values  between  48.2  mA  and  49.2 
mA.  The  laser  operation  of  /^,  =0.206  and  injection  current  1=47.87  mA  destroys  the  twin-mode 
locking  and  leads  to  a  new  single-mode  state  with  higher  output  power  6.279  mW.  For 
l>48.23mA  the  high  frequency  oscillation  occurs.  The  strongest  oscillation  shows  up  when 
1=48.39  mA  at  a  frequency  21.2  GHz  .  Further  increasing  injection  current  I  leads  the  weaker 
oscillation.  When  I>49.19niA,  the  stable  oscillation  disappears.  Thus,  the  injection  current  are 
functions  of  the  oscillation  frequency,  the  oscillation  amplitude  and  the  average  output  power 


(a) 
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T - 1 - 1 - 1 - r 


48.2  48.3  48.4  48.5  48.6  48.7  48.8  48.9  49  49.1  49.2 

injection  current  I  (mA) 


(c) 

Fig.  9  The  parameters  of  the  laser  operation  depend  on  the  injection  current  (I)  at  an 
external  cavity  length  Lext=5.25  mm  and  the  feedback  parameter =0.206  .  Data 
asterisks  are  the  simulation  data.(a)  Oscillation  frequency  versus  the  injection  current  (I). 

(b)  Average  output  power  versus  the  injection  current  (I),  (c)  Oscillation  Amplitude 

versus  the  injection  current  (I). 

According  to  the  simulated  results,  the  SEC  LD  operation  do  exist  some  instability 
regimes.  The  feedback  parameters  of  the  onset  and  offset  of  the  instability  depends  on  the  output 
power.  The  oscillation  frequencies  also  depends  on  the  output  power.  But  the  theoretical  critical 
feedback  value  in  equation  (27)  cannot  show  any  information  about  the  relation  between  and 
Power.  The  theoretial  frequency  value  in  equation  (28)  also  cannot  give  any  information  about 
the  relation  between  frequency  and  Power.  Those  should  be  investigated  in  more  detail. 

An  interesting  question  is:  what  is  the  highest  oscillation  frequency?  By  inserting 
k=C,Ito  into  equation  (24)  and  equation  (24)  ;r+  {zJq  =  -Arrsin^zio  into  equation  (29),  we  get 

v=  [-sin^zio]— •—  (42) 

where 

is  a  coupling  strength  coefficient  between  the  EC  and  the  LD.  From  equation  (42)  ,  we  can  see 
that  higher  and  shorter  LD  (means  shorter  r^)  lead  to  higher  oscillation  frequency. 

To  obtain  the  highest  oscillation  frequency,  we  also  need  to  set  -  sin  <I>q=\  in  equation 
(42).  Then  we  get 
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1 


assuming  t]=\.  Thus,  the  highest  frequency  is  equal  to  the  mode  spacing,  v=\ !  =  c !  2Ln^  , 

for  the  LD  without  the  EC. 

Whether  this  high  frequency  limit  can  really  be  achieved  and  where  the  stability  of  the 
high  frequency  regimes  really  are  still  need  to  be  investigated  in  more  detail. 

3.2  VCSEL  with  an  External  Cavity 


VCSEL  External  Cavity 


R1 


R2 


1 - + - 1 

T  Lcxt 

The  system  is  schematically  represented  in  Fig.  10  ,  wnere  R1  =0.999  (for  the  top  mirror), 
R2=0.995  (for  the  b(  Schematic  of  a  VCSEL  cavity  with  an  external  cavity  im.  The  other 
parameters  used  in  tae  louowmg  saiiuiauoii  are  siiown  in  i  aoie  in  . 


According  to  our  simulation,  the  result  of  simulation  agrees  with  analysis  only  when  the 
average  output  power  of  VCSEL  is  around  15.7  mW  *.  This  power  is  too  high  for  the  VCSEL 
opearion.  Hence  we  will  show  the  simulated  results  corresponding  to  the  output  power 
P=2mW/the  bottom  facet.  In  Fig.ll,  we  show  the  locations  of  strongest  high  frequency 
oscillations  in  the  parameter  space  consisting  of  the  feedback  parameter  (/^^,).and  the  length  of 
external  cavity  (Text).  The  graphical  representation  of  the  high  frequency  oscillations  (f)  as  a 
function  of  the  length  of  the  external  cavity  (Text)  is  shown  in  Fig.  12.  As  before,  the  dash  lines 
are  fits  to  simulated  results  marked  with  asterisks. 


TABLE  III 

List  of  Parameters  for  VCSEL 


Parameter 

Symbol 

Value 

Units 

scattering  loss  in  the  active  region 

3.6x10^ 

m"' 

carrier  number  at  transparency 

1.41372x10’ 

23 


output  power 

P 

2.653x10'*  X  5 

mW  /  the  bottom 

facet 

index  of  active  region 

3.5 

diameter  of  VCSEL’s 

r 

15 

um 

differential  gain 

Sn 

2.2x10* 

5'' 

linewidth  enhancement  factor 

a 

4 

laser  wavelength 

X 

966 

nm 

spontaneous  emission  lifetime 

2.0x10'® 

s 

photon  lifetime 

1.766x10''* 

s 

gain  saturation  coefficient 

K 

4.8 

w-' 

Fig.  1 1  A  VCSEL  with  a  SEC  operation  (Lext  ,fext).  The  lower  and  upper  solid  lines 
show  the  critical  value  at  and  group.  The  lower  and  upper  dash  lines  show 

the  strongest  oscillation  at  and  group  Curved  marked  with  represnet  the 
simulated  data.  The  parameters  for  the  strongest  oscillation  at  an  average  output  power 
around  15.7mW  .coincide  with  the  critical  parameters  for  the  onset  of  oscillation 
according  to  analysis.  The  range  of  parameters  for  oscillations  is  relatively  small. 

As  we  can  see,  shorter  laser  cavity  length  as  for  the  case  of  a  VCSEL  also  lend  to  twin¬ 
mode  locking  and  high  frequency  oscillations.  The  lower  dash  line  in  Fig.l  1  and  Fig.  12  show  the 

strongest  oscillation  just  before  switching  to  a  new  single-mode  state  for  and  the  upper 
ones  for  .  Increasing  /„,  causes  the  instability  to  switch  to  a  new  single-mode  state.  This 
process  repeats  if  f^,  continues  to  increase  (for  ,  N=0,l  ,2,3 . ) . 
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Next,  we  will  explain  the  instability  regime  of  VCSEL  operation  for  an  external  cavity 
length  Lext  =  3.75  mm  (the  instability  group  for  fj  only).  The  onset  of  instability  is 
fext  with  the  weak  intensity  of  sidemode.  Increasing  f^,  results  in  stronger  intensity  of 
sidemode  thus  stronger  oscillation.  When  /„,  reaches  0.071  as  show  in  Fig.  13  the  modulation 
depth  reaches  unity  (the  strongest  oscillation)  with  a  frequency  of  18ghz.  Increasing  further 
destroys  the  twin-mode  locking  instability  and  lead  to  a  higher  power  (~4.496mW)  in  a  single¬ 
mode  state.  These  changing  in  the  strength  of  the  oscillation  repeat  for  ,  N=0,l,2,3 . 
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Fig.  12  A  VCSEL  with  a  SEC  operation  (Lext  ,f).  The  lower  and  upper  solid  lines  show 
the  critical  beat  frequencies  of  the  and group.  The  lower  and  upper  dash 

lines  show  the  beat  frequencies  of  the  and  group  according  to  the  strongest 
oscillation  in  Fig.l  1 .  Curved  marked  with  represent  the  simulated  data. 
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Fig.  13  The  strongest  oscillation  in  Lext=3.75mm  and  the/^^,“  group.  The  feedback 
parameter  is  0.071.  The  oscillation  frequency  is  18  GHz. 


25 


Next,  we  will  show  the  dependences  of  the  oscillation  frequency  (f),  average  power  (P) 
and  oscillation  amplitude  on  the  injection  current  (I=J/e).  The  external  cavity  length  is  fixed  at 
Lext=3.75mm  and  the  feedback  parameter  is  fixed  at  /„,  =0.071.  The  other  parameters  are 
shown  in  Table  III.  In  Table  IV  and  Fig.  14  (a),(b),(c)  ,  we  can  see  the  variations  of  these 
quantities  for  injection  current  values  between  138.8mA  and  193.8mA.  The  laser  operation  for 
=0.071  and  injection  current  1=1 38.83mA  detroys  the  twin-mode  locking  and  leads  to  a  new 
single-mode  state  with  higher  output  power  4.329m W.  For  I>1 3 8.99mA  the  high  frequency 
oscillation  occurs.  The  strongest  oscillation  shows  up  when  1=1 39.07mA  at  a  frequency  18 
GHz.  Further  increasing 

Table  IV 

List  of  frequency,  average  power  and  oscillation  amplitude  on  the  change  of  the  injection  current. 


injcaion  curi«nt  I  (mA) 


(a) 


26 


2.5 


2.4 


1.9 


1.8' - ' - ' - ' - ■ - - - - - - - ■ - ‘  . ^ 

138.8  138.9  139  139.1  139.2  139.3  139.4  139-5  139.6  139.7  139.8 

injection  cufxent  1  (mA) 

(b) 


1.5 

I 

I  > 
J 

i*- 


138.8  138.9  139  139.1  139-2  139-3  139-4  139-5  139-6  139-7  139-8 

injection  cutrent  I  (mA) 

(c) 

Fig.  14  The  parameters  of  the  laser  operation  depend  on  the  injection  current  (1)  at  an 
external  cavity  length  Lext=3.75mm  and  the  feedback  parameter  =0.071  .  Data 
marked  asterisks  are  the  simulation  data,  (a)  Oscillation  frequency  versus  the  injection 
current  (I),  (b)  Average  output  power  versus  the  injection  current  (I),  (c)  Oscillation 
Amplitude  the  injection  current. 

injection  current  I  leads  the  weaker  oscillation.  When  I>  139.63mA,  the  stable  oscillation 
disappears. 

Due  to  the  simulated  results  in  Fabry-Perot  lasers  and  VCSEL,  we  can  conclude  that:  the 
SEC  LD  exhibits  some  regimes  of  high-frequency  self-oscillations  in  some  range  of  feedback 
parameters.  For  fixed  Lext  and  ,  the  change  of  the  injection  current  will  cause  the  change  of 
the  oscillation  frequency,  the  oscillation  amplitude  and  the  average  output  power. 
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4.0  Conclusion 


We  have  presented  analytical  and  numerical  investigation  of  the  instability  regimes  of  a 
SEC  LD  .  The  instability  is  caused  by  self  twin-mode  locking  of  a  short  external  cavity  (a  few 
millimeters)  laser  diode.  As  we  have  proved  theoretically,  these  high-frequency  oscillations 
correspond  to  an  equalization  of  the  cavity  loss  for  the  main  mode  (the  lasing  frequency)  and  the 
sidemode.  The  phenomena  of  high-frequency  oscillation  is  caused  by  a  high  level  of  coupling 
between  the  laser  diode  and  the  external  cavity  and  thus  are  compound-laser-cavity  effects. 

The  key  results  are  as  follows:  First,  a  laser  diode  operation  with  a  short  external  cavity  (a 
few  millimeters)  in  some  ranges  of  feedback  parameters  exist  some  regimes  of  high-frequency 

self  twin-mode  oscillations.  Second,  in  high  frequency  oscillation  of  every  (N=0,l,2 . ) 

regime,  increasing  leads  to  higher  modulation  index  and  higher  oscillation  frequency.  For 
certain  the  strongest  oscillation  amplitude  shows  up  (the  modulation  depth  reaches  unity). 
Increasing  destroys  the  twin-mode  locking  instability  and  lead  to  a  higher  power  single¬ 
mode  state.  These  instability  repeat  for  every  /^,  ^  regime. 

In  the  simulation  results,  changing  injection  number  (J)  (or  injection  current  (I))  will 
change  the  average  output  power.  In  addition  to  changing  the  output  power,  it  also  changes  the 
critical  feedback  parameter  which  required  for  high-frequency  oscillations  and  the  oscillation 
frequency.  But  in  equations  (26)  and  (29),  the  critical  feedback  parameters  and  the  beat 
frequencies  do  not  depend  on  either  the  injection  number  or  the  output  power.  The  validity  of 
these  questions  should  be  investigated  further. 

In  numerical  simulation,  we  used  two  kinds  of  laser  cavities,  Fabry-Perot  laser  and 
VCSEL.  The  numerical  and  analytical  results  demonstrate  that  a  SEC  LD  operation  do  exhibit 
high-frequency  self  twin-mode  locking  oscillations  in  some  ranges  of  feedback  parameters. 
These  should  be  proved  by  experiments. 

This  laser  operation  of  a  SEC  LD  can  be  used  as  tunable  microwave  oscillators  with 
frequency  modulation  capacities.  These  device  can  be  very  useful  for  all-optical  high-speed 
microwave  fiber-optic  transmitters. 
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Appendix  A 


To  find  the  relation  between  output  power  (P)  and  photon  numbers  (S),  we  should  discuss 
the  two  main  sources  of  loss  in  optical  resonators  first.  The  two  most  common  loss  mechanisms 
in  an  optical  resonator  are  as  follows: 

1.  Losses  resulting  from  absorption  and  scattering  in  the  laser  medium  between  the  mirrors,  a, 
is  the  absorption  coefficient  of  the  medium  corresponds  to  photons  which  are  lost  within  the 
laser  cavity. 

2.  Losses  arising  from  nonperfect  reflections  at  the  mirrors.  Reflection  loss  is  unavoidable, 
since  without  some  transmission  no  power  output  is  possible.  Even  we  design  the  mirrors 
with  the  highest  possible  reflectivities,  some  residual  absorption  and  scattering  reduce  the 
reflectivity  to  somewhat  less  than  100  percent. 

3. 

The  threshold  gain  denoted  g  =  g,*  ,  equals  to  an  effective  overall  distributed-loss 
coefficient, 

g,h  =^s+ 

where  the  quantity 


1  w  1 

a„,  = - ln( - 

“  2L 


(A-2) 


represents  a  fictitious  equivalent  loss  per  unit  length  due  to  useful  escape  of  light  through  one  of 
the  two  cavity  mirrors.  The  reflectances  of  the  two  mirrors  are  R1  and  R2. 


The  external  quantum  efficiency  of  the  LD,  rj^,,  is  the  fraction  of  injected  electrons  that 
get  converted  to  usefiil  output  photons.  The  quantity  77,  is  the  internal  quantum  efficiency,  the 
fraction  of  injected  electrons  that  produce  stimulated  photons,  not  all  of  which  will  reach  the 
output  of  the  LD.  The  ratio  between  the  external  and  the  internal  quantum  efficiency  is  given  as 


'Hexl  _  _ ^1-^2  _ 

Vi  g,h  g,h  +  < 


(A-3) 


Therefore,  total  emitted  power  due  to  stimulated  emission  is 


=  (/2  v)  •  •(  #  of  photons  due  to  stimulated  emission  per  unit  time)  (A-4) 

Vi 


that  is. 
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Ik 

=  ihv)' 


li 


Ve.,  1 


■s 


(A-5) 

(A-6) 


Vi  •z'pA 

where  •  S  is  the  photons  generated  per  unit  time  by  stimulated  emission  and 


R. 


1 


<'ph 


(A-7) 


is  above  the  threshold.  The  photon’s  lifetime,  ,  is  defined  as 


1  /  C  X  /  ^ 


(A-8) 


where  is  the  effective  group  refractive  index  and  hv  \s  the  energy  of  each  photon. 
From  equations  (A-3)  and  (A-8),  the  total  emitted  power  will  be 


P,o,=ihv)^ 


oc. 


n„ 


Then, 


P,oi  =ihv)-a„,-  —  S 


(A-9) 


(A-10) 


With  the  expression  for  given  by  equation  (A-2),  one  finally  gets 


rig  2  L  y?,  R2 


(A-11) 


where  P,^,  is  the  total  power  emitted  from  both  facets  (or  both  mirrors).  Equation  (A-1 1)  is  the 
relation  between  power  (P)  and  photon  numbers  (S). 

For  a  symmetric  laser  diode,  a  total  output  power  is  equal  to  twice  of  the  output 
power  per  facet  (P  /  per  facet  =  1/2  i^„,) . 
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Appendix  B 

Small  signal  analysis  of  the  EC  LD  far  above  threshold  can  be  performed  by  linearizing 
(1)  and  (2)  vv^ith  respect  to  small  derivations  around  the  stationary  solutions  and  neglecting  the 
noise  terms.  By  inserting  S  =  jco  into  equation  (32)  ,  we  get 

D(o))  =  io/  -ty^[y  +  2A:(l-e'‘"^)cos^] 

+/c[<u^^(l -e"“'’)(cos^^- orsin^)]  (B-1) 


with  e(  r)  =  1  -  (j)  =  Qr  .  Then, 

D{co)  =  ico^  -  (o^^y  +  lKil  -  cos cot+  i sin 6>r]cos  ^ 

-ico^o)^  +  Ar^(l-2(cos6Jr+/sin0r)  +  cos2fyr+/sin2n>r)] 

(1  -  cos  (Dt- i  sin  (ot)(cos^-  asin  (j))^  (B-2) 

For  having  a  zero  of  D{(d)  on  the  imaginary  axis  (the  boundary  of  the  stability),  we  set 
the  real  part  and  imaginary  part  of  D{(o)  equal  to  zero  (Re[Z)((»)]=0  and  Im[D(a))]=0)  .  The 
real  part  of  D{co)  gives 

Re[£)((ij)]=  -6)^  [y  +  2/c(l  -  cos  cor)  cos  (-2/  sin  cot+  i  sin  20^7)] 

+Kco^  (1  -  cos  coT)(cos(p-  (zsin  (fi) 

=  [y  +  2a:(1  -  cos  cot)  cos  (f\  -  co[k^{-2  sin  cot+  2  sin  cor-  cos  cor)^ 

+Ka)^  (1  -  cos  fyT)(cos  (p-  asin 

=  -o)^Y-2(o^k{\-  cos  (OT)cos<j)+  coK^  (-2  sin  cot+  2  sin  cor- cos  cor) 

+ K(o^{\- cos  (OT){Q,os(p- asm  (f)) 

For  Re[D(iu)]=0,  we  get 

0=  -nj^y-2nj^A:(l  -cos6Jr)cos^+  a)A^(-2sina;r+2sin<yT-cosajr) 

+/C£y/(l-cos(yr)(cos^-asin^)  (B-3) 


Dividing  equation  (B-3)  by  {kco^)  ,  we  obtain 

.  CO^  ^  0)^  ^  COK  . 

0  = - 2 —  2 - j  ■  /c(l  -  cos  cot)  cos  f-2 - Y 

r 

+  (cos  a  sin  (j)){\  -  cos  cot) 

Then, 


(B-4) 
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(1  -  2  •  -^)(1  -  cos  cot)  cos  ^  -  orsin  -  cos  cor) 


CO 


(B-5) 


„  COK  .  \  ^ 

=  2 - r  •  sm  car- (1  -  cos  oor)  + 


a. 


KCO„ 


Dividing  equation  (B-5)  by  (1  -  cos  cox)  gives 


(1-2 - y)  ^  -  asm  ^  =  2 - r-  •  sm  (ot+ 


CO, 


CO, 


KCOj  1  -  COS  COT 


(B-6) 


This  is  the  boundary  condition  from  the  real  part  of  D{co) 
The  imaginary  part  of  D{co)  gives  , 


Im[D(ro)]-  ico^  -co^[2iKS\ncoTCOS(fi\-ico[co^  +  k^{\-2coscot+cos2cot)] 

(B-7) 


+K\-co^i  sin  cot{cos^- asin  ^)] 
Setting  Im[D(6j)]=0  gets 


0)^  -  co^[2KsmcoTQ.os(l)\- co[co^  +  K^{\- 2  cos  COT +cos2 cot)) 
-h^co^  sm.o)T{cos<l)- asm  <f))]-0 


(B-8) 


Then, 


-  coco^  =  2k(o^  smcoTCos(f)+  cot^  {\- 2  cos  cor  +  cos2cot) 
+  Kco^  sin  a)T{cos^-  asin  (p) 


or 


co{co^  -  co^ )  =  2  Kco^  sin  cotcos  cj)  +  cok^  (1-2  cos  cot  +  cos 2  cot) 
+  Kco^  sin  cot{cos  (p  -  asin  cp) 


i.e. 


co{co^  -  co^)  =  coi^{\-2coscot+  cos2<yr) 


■\-KC0^  sincyr{  [l-2(— )^] cos  asin  (Zi) 


(B-9) 


CO, 


By  inserting  equation  (B-6)  into  RHS  of  equation  (B-9) ,  we  obtain 


co{co^  -  co^ )  =  coi^  (1-2  cos  COT+  cos2cot) 


COK 


+  KCOjf  sm  cot{2  — Y  ■  sm  cot+ 


co^ 


(B-10) 


CO. 


KCO/  1  -  COS  COT 
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Divided  by  o),  we  get 


ai^  -0}/  =  k^{\-2coso)t+cos2o}t) 

2  •  ^  •  o)  1  , 

+  KCOj^  sin  — -  •  sin  - ; - J 


CO, 


KO)j  1  -  COS  CDT 


Then, 


o/  -  0)/  =  ic^  {1-2  cos  COT+  cos2o)t) 
■2k^  ■ 


.  2  co  cOg  sin^ur 
+  1K^  -sm  COT-  ^  - 


COj  1  -  cos  COT 


CO, 


With  ooj  =  — —  ,  the  damping  frequency  of  relaxation  oscillation  of  LD  without  EC  ,  we  get 


co^  -  00^  =  /c^(l  -2cos6JT+cos2(i)r+2sin^  cot) 

^  .  COT  COT 

2  2  sin - cos  — 

_oo^ _ 2  2 

(^)  2sin^^ 

Y  ^ 


(B-11) 


Then, 


CO  -cOn  =  ^(2(l-cos(yr)  +  cos2tyr-(l-2sin  ^yr))  -  oy  •  cot — 


Using  cos 2617=  1 -2sin^  leads 

£9^ =4A:^sin^ -^-fiiy  cot-^  .  (B-12) 

Therefore,  equations  (B-6)  and  (B-12)  show  the  stability  boundary. 

The  state  becomes  unstable  when  a  zero’ of  D{co)  passes  the  imaginary  axis.  For  the 
stable  solutions  the  zeros  of  D{co)  should  lie  on  the  left  side  of  the  imaginary  axis.  Then  we  have 
the  stability  solutions  in  equation  (33)  and  (34)  . 
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Appendix  C 


/♦  C  PROGRAM  */ 

#mclude  <stdlib.h> 
#mclude  <stdio.h> 
#mclude  <math.h> 


#define 

#define 

#defme 

#defme 

#define 

#defme 

#defme 

#defme 

#defme 

#defme 

#defme 

#defme 

#defme 

#define 

#defme 

#defme 


KP 

R1 

R2 

L 

NG 

C 

alfs 

TSP 

NS 

GN 

A 

DL 

DH 

J 


4.8  /*  gain  saturation  coefficient  *! 

0.999  /*  facet  reflectivity  */ 

0.995  /*  facet  reflectivity  */ 

1  .Oe-6  /*  laser  cavity  length  *! 

3.5  /*  refractive  index  */ 

3.0e8  /*  light  velocity  */ 

3.6e3  /*  scattering  loss  coefficient  */ 

2.0e-9  /*  spontaneous  emission  lifetime  */ 

1 .4 1 372e9  /*  carrier  number  to  reach  zero  gain  */ 


2.2e3 
4.0 
0.0 
4e-8 
0.87el8 
INTERVAL  1000 
STEPS  4000000 


/*  gain  coefficient  */ 

/*  linewidth  enhancement  factor  */ 

/*  simulation  initial  time  *! 

/*  simulation  final  time  10ns  */ 

/*  injection  number  */ 

/*  data  printout  interval  */ 
/*  steps  number  in  (0,10ns)  */ 


void  odelongO; 
void  rk5(); 
void  diff(); 
void  FB(); 
void  nrerrorO; 
void  ffee_vector(); 
void  derivsO; 
double  *vector(); 

double  FEXT, TEXT, TP, alpha, P; 
int  delay; 


main() 

{ 

long  nvar,  i; 
double  xl,  x2,  hi,  *y; 


alpha  =  alfs-log(Rl*R2)/(2.0*L); 

TP  =  NG/(C*alpha); 

P=  -0.5*(6.63e-34*C/966e-9)*(log(Rl*R2)/(2.0*L))*(C/NG);  /♦wavelengh  966nm*/ 
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nvar=3; 

xl=DL; 

x2=DH; 


hl=(x2-xl)/STEPS;  /*  step  size  */ 

y=vector(l,  nvar); 
y[l]=1000.0;  /*  y[l]  1(0)  */ 

y[2]=0.0;  I*  y[2]  Q  */ 

y[3]=NS;  I*  y[3]  N(0)  */ 

printf("Input  the  value  of  FEXT  -=>"); 
scanf("%lf',&FEXT); 
printfC'Input  the  value  of  LEXT  -=>"); 
scanf("%lf',&LEXT); 

delay=2.0*LEXT/C/hl; 


printfC  P=%e  \n",  P); 

printf("  J=%e  \n",  J); 

printf("  delay=%d  \n",  delay); 

printfC  FEXT=%1.4f  \n",FEXT ); 

printfC  LEXT=%2.2fmm  \n",  LEXT*le3); 


odelong(y,  nvar,  xl,  x2,  hi,  derivs,  rk5); 


} 


void  nrerror(error_text)  /*  error  handle  */ 
char  error_text[]; 

{ 

void  exit(); 

fprintf(stderr,  "  numerical  run-time  error  ...  \n"); 
fprintf(stderr,  "%s  \n",  error_text); 
fprintf(stderr, "  ...exiting  to  system  ...  \n");  exit(l); 

} 

double  *vector(nl,  nh) 
long  nl,  nh; 

/*  allocates  a  double  vector  with  range  [  1  ..  nvar]  */ 

{ 

double  *v; 

v=(double  *)malloc((unsigned)  (nh-nl+l)*sizeof(double)); 
if(!v) 

nrerror("allocation  failure  in  vector()"); 
return  v-nl; 

} 
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void  free_vector(v,  nl,  nh) 
double  *v; 
long  nl,  nh; 

/*  frees  a  double  vector  allocated  by  vector()  */ 

{ 

free((char*)  (v+nl)); 

} 


void  derivs(x,  y,  v,  dydx) 

I*  yl=S,  y2=Q  ,  y3=N 
differential  equation  */ 

double  X,  y[],  v[],  dydx[]; 

{ 

double  CL,  TO,  K,  tempi,  ss,  Q; 

. /*  double  CL,  TO,  K,  tempi ,  ss,  Q,  R1 ;  */ 

CL=0.5*(l-R2)/sqrt(R2); 

TO=2*L*NG/C; 

K=2*CL*sqrt(FEXT)/TO; 

templ=GN*(y[3]-NS); 

ss=temp  1  *  ( 1  -ICP*P*y[  1  ])-  1/TP; 

Q=  -atan(A)+y[2]-v[2]; 


dydx[  1  ]=ss*y  [  1  ]+2*K*  sqrt(y  [  1  ])*sqrt(v[  1  ])*  cos(Q); 
dydx[2]=0.5  *  A*temp  1  -K*sqrt(v[  1  ])*sin(Q)/sqrt(y[  1  ]); 
dydx[3]=J-y  [3]/TSP-temp  1  *  ( 1  -KP*P*y  [  1  ])*y  [  1  ]; 

} 


void  odelong(ystart,  nvar,  xl,  x2,  hi,  derivs,  rk5) 
double  ystart[],  xl,  x2,  hi; 
long  nvar; 

void(*derivs)();  /*  ANSI:  void(*derivs)(double,double  ‘'',double  *);  */ 
void(*rk5)(); 

/*  runge-kutta  integrate  starting  values  ystart[l,  nvar]  from  xl  to  x2, 
hi  should  be  set  as  stepsize. 

derivs  routine  for  calculate  the  right-hand  side  derivative,  while  rk5 
is  the  name  of  the  stepper  routine  to  be  used  *! 


{ 

longi,  k,  step; 
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double  X,  h; 

double  *y,  *v,  *dydx,  *yout,  *el,  *e2,  *vector(); 
void  nrerrorO,  free_vector(); 

y=vector(l,  nvar); 
v=vector(l,  nvar-1); 
dydx=vector(  1 ,  nvar); 
y  out=vector(  1  ,nvar) ; 
h=hl; 


el=vector(l,  STEPS+1);  /*  data  storage  for  S(t)  *! 
e2=vector(l,  STEPS+1);  /*  data  storage  for  phase  */ 


printfC  h=%5.10e\n",  h); 
printfC  max-step=%d  \n",  STEPS); 
printfC  step-delay=%d  \n",  delay); 


for  (i=l;  i<=nvar;  i++) 
y[i]=ystart[i]; 
x=xl; 
v[l]=0.0; 
v[2]=0.0; 

for  (step=l;  step  <=  STEPS;  step++)  { 

el[step]=y[l]; 
e2[step]=y[2]; 
if(step  <=  delay ) 

{ 

v[I]=0.0; 

v[2]=0.0; 

} 

else 

{ 

v[l]=el  [step-delay]; 
v[2]=e2  [step-delay] ; 

} 

(*derivs)(x,  y,  v,  dydx); 

rk5(y,  v,  dydx,  nvar,  x,  h,  el,  e2,  step,  yout,  derivs,  diff,  FB); 
x=x+h; 

for  (i=l;  i<=nvar;  i++) 
y[i]=yout[i]; 


} 


el[step]=y[l]; 

e2[step]=y[2]; 
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x=xl; 


for  (i=0;  i<STEPS/INTERVAL;  i++)  /*  print  out  the  datas  */ 

{ 

k=INTERVAL*i+l; 

I*  printfC  x=%e  \t",  x);  */ 
printf("  %e  \t",  x); 
printfC  %5.3f\n",  el[k]*P*le3); 
x=x+INTERVAL*h; 

} 

/*  printfC  x=%e  \t",  x);  */ 

printfC  %e  \t",  x); 

printfC  %5.3f\n",  el[k]*P*le3); 

I*  fprmtf(fp,"  %e  \t",  x);  */ 

C  fi)rintf(fi),"  %5.3f\n",  el[k]*3.6e-5);  */ 

ffee_vector(el,  1,  STEPS+1); 
ffee_vector(e2,  1,  STEPS+1); 
free_vector(y,  1 ,  nvar); 
free_vector(v,  1,  nvar-1); 
free_vector(yout,  1,  nvar); 
ffee_vector(dydx,  1  ,nvar); 


} 


void  rk5(y,  v,  dydx,  n,  x,  h,  el,  e2,  step,  yout,  derivs,  diff,  FB) 
double  y[],  v[],  dydx[],  x,  h,  yout[],  el  [],  e2[]; 

void  (*derivs)();  /*  ANSI:  void  (*derivs)(double,double  *,double  *);  */ 
void  (*diff)(); 
void  (*FB)(); 
long  n,  step; 

/*  given  values  for  n  variables  y[l,n]  and  their  derivatives  dydx[l,n]  known 
at  X,  use  the  fifth  order  Runge-Kutta  method  to  advance  the  solution  over 
an  interval  h  and  return  the  incremented  variables  as  yout[l,n].  which  need 
not  be  a  distinct  array  from  y.  the  user  supplies  the  routine  derivs(x,y,dydx) 
which  return  derivatives  dydx  at  x.  */ 

{ 

long  i; 

double  *vl,*Fl,  *F2,  *F3,  *F4,  *F5,  *F6,  *dydxt,  *ytemp,  *vector(); 
double  *S,*T; 
void  free_vector(); 

Fl=vector(l,  n); 

F2=vector(l,  n); 

F3=vector(l,  n); 

F4=vector(l,  n); 

F5=vector(l,  n); 
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F6=vector(l,  n); 
ytemp=vector(l,  n); 
dydxt=vector(l,  n); 

S=vector(l,5); 

T=vector(l,5); 

vl=vector(l,l); 

if  (step>delay) 

{ 

(*diffXh,step,el,S); 

(*diff)(h,step,e2,T); 

} 


for  (i=l;  i<=n;  i++) 

Fl[i]=h*dydx[i];  /*  first  step  */ 
for(i=l;  i<=n;  i++) 
ytemp[i]=y[i]+0.5*Fl[i]; 


if  (step>delay) 

{ 

(*FB)(0.5*h,step,S,el,vl); 

v[l]=vl[l]; 

(*FB)(0.5*h,step,T,e2,vl); 

v[2]=vl[l]; 

} 


(♦derivs)(x+0.5*h,  ytemp,  v,  dydxt);  /*  second  step  */ 


for(i=l;  i<=n;  i++) 

F2[i]=h*dydxt[i]; 
for(i=l;  i<=n;  i++) 

ytemp[i]=y[i]+0.25*Fl[i]+0.25*F2[i]; 
if  (step>delay) 

{ 

(♦FBX0.5*h,step,S,el,vl); 

v[l]=vl[l]; 

(*FB)(0.5*h,step,T,e2,vl); 

v[2]=vl[l]; 

} 

(*derivs)(x+0.5*h,  ytemp,  v,  dydxt); 


for(i=l;  i<=n;  i++) 

F3[i]=h*dydxt[i];  /*  third  step  */ 

for(i=l;  i<=n;  i++) 

ytemp[i]=y[i]-F2[i]+2.0*F3[i]; 
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if  (step>delay) 

{ 

v[  1  ]=e  1  [step-delay+ 1  ] ; 
v[2]=e2[step-delay+l]; 

} 

(*derivs)(x+h,  ytemp,  v,  dydxt); 
for  (i=l;  i<=n;  i++) 

F4[i]=h*dydxt[i];  /*  fourth  step  */ 

for(i=l;  i<=n;  i++) 

ytemp[i]=y[i]+7.0*Fl[i]/27.0+10.0*F2[i]/27.0+F4[i]/27.0; 

if  (step>delay) 

{ 

(*FB)(2.0*h/3.0,step,S,el,vl); 

v[l]=vl[l]; 

(♦FB)(2.0*h/3.0,step,T,e2,vl); 

v[2]=vl[l]; 

} 

(*derivs)(x+2.0*h/3.0,  ytemp,  v,  dydxt); 
for  (i=l;  i<=n;  i++) 

F5[i]=h*dydxt[i];  /*  fifth  step  */ 

for(i=l;  i<=n;  i++) 

ytemp[i]=y[i]+28.0*Fl[i]/625.0-0.2*F2[i]+546.0*F3[i]/625.0+54.0*F4[i]/625.0-378.0*F5[i]/625.0; 

if  (step>delay) 

{ 

(♦FB)(0.2*h,step,S,el,vl); 

v[l]=vl[l]; 

(*FB)(0.2*h,step,T,e2,vl); 

v[2]=vl[l]; 

} 


(*derivs)(x+0.2*h,  ytemp,  v,  dydxt); 
for(i=l;  i  <=n;  i++) 

F6[i]=h*dydxt[i];  /*  sixth  step  */ 

for  (i=l;  i<=n;  i++) 

yout[i]=y[i]+Fl[i]/24.0+5.0*F4[i]/48.0+27.0*F5[i]/56.0+125.0*F6[i]/336.0; 

free_vector(dydxt,  l,n); 
free_vector(ytemp,  1,  n); 
free_vector(F6,  1,  n); 
free_vector(F5,  1,  n); 
free_vector(F4,  1 ,  n); 
free_vector(F3,  1,  n); 
free_vector(F2,  1,  n); 
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free_vector(Fl,  1,  n); 

free_vector(S,  1 , 5); 
free_vector(T,  1 , 5); 
free_vector(v  1,1,1); 

} 


void  FB(h,step,Sl,S2,vl) 
double  h,Sl[],S2[],vl[]; 
long  step; 

/*  Taylor-series  expansions 
Found  the  feedback  value  between  t  and  t+h  */ 

{ 


vl[l]=S2[step-delay]+Sl[l]*h+Sl[2]*h*h/2.0+Sl[3]*h*h*h/6.0+Sl[4]*h*h*h*h/24.0+Sl[5]*h*h*h*h*h/120.0; 


} 


void  diff(h,step,S3,S4) 
double  h,S3[],S4[]; 
long  step; 


/*  numerical  differentiation  : 

if  step<=delay+3  :  using  Forward-Difference  Expressions  with  Error  of 
Order  h*h 


if  step>delay+3  :  using  Central-difference  Expressions  with  Error  of 
Order  h*h*h*h 


*! 


{ 

long  i; 

i=step-delay; 

if(step>delay+3) 

{ 

S4[l]=(-S3[i+2]+8*S3[i+l]-8*S3[i-l]+S3[i-2])/(12*h); 

S4[2]=(-S3[i+2]+16*S3[i+l]-30*S3[i]+16*S3[i-l]-S3[i-2])/(12*h*h); 

S4[3]=(-S3[i+3]+8*S3[i+2]-13*S3[i+l]+13*S3[i-l]-8*S3[i-2]+S3[i-3]) 

/(8*h*h*h); 

S4[4]=(-S3[i+3]+12*S3[i+2]-39*S3[i+l]+56*S3[i]-39*S3[i-l]+12*S3[i-2] 

-S3[i-3])/(6*h*h*h*h); 
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else 


{ 

S4[l]=(-S3[i+2]+4*S3[i+l]-3*S3[i])/(2*h); 

S4[2]=(-S3  [i+3]+4*  S3  [i+2]-5*  S3  [i+ 1  ]+2*  S3  [i])/(h*h)  ; 

S4[3H-3*S3[i+4]+14*S3[i+3]-24*S3[i+2]+18*S3[i+l]-5*S3[i])/(2*h*h*h); 

S4[4]=(-2*S3[i+5]+ll*S3[i+4]-24»S3[i+3]+26*S3[i+2]-14»S3[i+l]+3*S3[i]) 
/(h*h*h*h) ; 

} 


42 


References 


1.  BJARNE  TROMBORG  and  JESPER  MORK,  “Nonlinear  injection  locking  dynamics  and 
the  onset  of  coherence  collapse  in  external  cavity  lasers,”  IEEE  J.  Quantum  Electron.  , 
vol.  QE-26,  pp. 642-654,  Apr.  1990 

2.  J.  Mork  ,  B.  Tromborg  ,  and  J.  Mark  ,  “  Chaos  in  semiconductor  lasers  with  optical 
feedback  :  Theory  and  experiments,”  IEEE  J.  Quantum  Electron.  ,  vol.  28,  pp.  93-108,  Jan. 
1992. 

3.  R.  Lang,  and  K.  Kobayashi,  “External  optical  feedback  effects  on  semiconductor  injection 
laser  properties,”  IEEE  J.  Quantum  Electron,  vol.  16  ,  pp.347-355.  Mar.  1980. 

4.  Alexei  A.  Tager  and  Klaus  Petermann,  “High  frequency  oscillations  and  self-mode  locking  in 
short  external-cavity  laser  diodes,”  IEEE  J.  Quantum  Electron. ,  vol. 30  ,  pp.1553-1561,  July 
1994. 

5.  Shtengel,  G.,  Temkin.  H.,  and  Swirhun,  S.E.  ,  “High-spped  vertical-cavity  surface  emitting 
laser,”  IEEE  Photonics  Technol.  Lett. ,  1993,  5,  (12),  pp.  1359-1362 

6.  Petermann ,  Laser  diode  modulation  and  noise  .  Kluwer  Academic  Publishers,  1991 . 

7.  Tromborg,  J.  H.  Osmundsen,  and  H.  Olsen,  “Satbility  analysis  for  a  semiconductor  laser  in 
an  external  cavity,”  IEEE  J.  Quantum  Electron.,  vol.  QE-22,  pp.294-301,  Feb.  1986  . 

8.  H.  Henry  and  R.  F.  Kazarinov,  “Instability  of  semiconductor  lasers  due  to  optical  feedback 
from  distant  reflectors  ,”  IEEE  J.  Quantum  Electron. ,  vol.  QE-22  ,  pp.  294-301,  Feb.  1986. 

9.  A.  Tager  and  B.  B.  Elenkrig,  “Stability  regimes  and  high-frequency  modulation  of  laser 
diodes  with  short  external  cavity,”  IEEE  j.  Quantum  Electron,  vol.  QE-29  ,  pp. 2886-2890, 
Dec.  1993  . 

10.  William  H.  Press,  Brian  P.  Flannery,  Saul  A.  Teukolsky  and  William  T.  Vetterling, 
Numerical  Recipes  :  The  art  of  scientific  computing  .  Cambridge  University  Press  1986  . 

11.  David  Kincaid  and  Ward  Cheney,  Numerical  Analysis  .  Brooks/Cole  Publishing  company. 
Pacific  Grove,  California,  1991  . 


43 


12.  Schunk  and  K.  Petermann,  “Numerical  analysis  of  the  feedback  regimes  for  a  single-  mode 
semiconductor  laser  with  external  feedback,”  IEEE  J.  Quantum  Electron.,  vol.  24,  pp.  1242- 
1247,  July  1988. 


44 


OU.S.  GOVERNMENT  PRINTING  OFFICE:  1997-509-127-47134 


MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 

a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Material 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence, 
reliability  science,  electro-magnetic  technology,  photonics,  signal 
processing,  and  computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


